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Abstract 


Properties  of  matrices  of  the  form  H=^(a.-  b<)~  j 


are  investigated.   If  G  =  H~  =   c,  .   then  explicit 


i  u3 

formulas  are  given  for  1°  5/) »  and  the  row  and  column 
sums  of  G.   These  extend  formulas  previously  given  by 
Collar  and  Smith  for  the  generalized  Hilbert  matrix. 
For  H  symmetric  and  positive  definite  the  smallest 
eigenvalue  of  H  and  its   P-  condition  are  estimated. 
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NYO-2539 


On  the    Inversion   of  Certain  Matrices 


1.      Introduction. 

Let      a,,ap,...,a     ;     b,,b2,...,b        be     2n     distinct, 
but   otherwise   arbitrary,    complex  numbers.      From  the   matrix 
H,    of   order     n, 

(1.1)  H  =     — — -^ — ;;        with      1   <   i,      j   <    n 

and  let   G  =  H~   =  )c..j,   (The  indices   i,j,k  will  range 
from   1  to  n  unless  it  is  specified  otherwise.)   If,  for 
some  constant   p, 

(1.2)  a.-b.  =  i  +  j-l+p^0 

u 

then  H  is  a  segment  of  the  well  known  generalized  Hilbert 
matrix,  and  in  this  case  formulas  for  c..  have  been  given 
by  Savage  and  Lukacs  [l\.] ,   Smith  [6]  and  Collar  [l].   For 

(1.3)  a±   -  b.    =     i  -  j  +  p 

Linfoot  and  Shepherd  [31  and  Collar  [2]  give  formulas  for 

c. .   and  in  both  cases  Collar  exhibits  diagonal  matrices 

1J  m 

T 

D,  K  such  that   G  =  DH  K.   Collar  [2]  and  Smith  also 

evaluate  the  quantities  \         u.,  .,  Y~      c.  .. 

1,7       j 


-  k   ~ 


These  authors  make  use,  in  most  cases,  of  the  formula 
for  the  determinant  [£] 


TT(a1-ak)(bk-  b.) 

(1.1;)  det  H  =  -^ 

THa,  -  bk) 

and/or  require  the  evaluation  of  certain  involved  series,, 

The  formulas  of  Collar  and 'Smith  are  extended  here  to 
the  general  case  (1.1).   The  method  to  be  used  does  not 
depend  on  (1.1;)  but  simply  on  Lagrange's  interpolation 
formula .   Indeed  (1.1;)  comes  out  as  a  by-product  of  formula 
(3.3)  given  below. 

2.   Formula  for  the  inverse. 


Let 

(2.1)      A(x)  =  TTU-a.)  ,    B(x)  =T|~U-b.) 

i      x  i      x 

and  denote  the  fundamental  polynomials  of  the  Lagrangian 
interpolation  corresponding  to  the   a.   and  b.,  respectively, 
by 

where  prime  denotes  differentiation.  We  then  have 
Theorem  1.   The  elements  of  G  are  given  by 

(2.3)  c^  =  (aj-bi)Aj(bi)B1(aJ.) 
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and  if  H   is  symmetric 

(24)     c±,  =  (aj  -  b1)AJ(bi)A1(bj). 

Proof?  For  any  polynomial  p(x)   of  degree  n-1   we 
may  write 

(2.5)  pU)  =  EZ  p(»4)  Mx> 

i  x        x 


or 

(2, 

.6) 

where 

(2, 

►7) 

m-^^ 


a. 


P(aj) 
Ci     ~     3157) 

Pk(&i) 
Now   let      p(x)    ■  -^(x)A(bk)    =  pk(x)      and     ckl  =  A>  (a    )    . 

Then  from      (2.6)   we   obtain 


(2.8)  Vx)A(V     =     5—    _^ki 

A(xJ  ^j-     aj_- 


x   * 


If  we  set     x  =  b,      then,   since     ^(b^)   =  Bw,   the  Kronecker 
delta,   we   get    that 

Thus      cw      gives   the  desired   inverse  element,    that   is, 


'ki 


B.  (a,)A(bJ 

^ki  =  -        A'(a4)  k-Ca1-^BkU1)A1(V 


—      o     ™* 


If  H   is  symmetries   a.-b.  =  a.-b.,  it  follows  that 
B.  (a  . )  =  Aj(b  . )   and  the  theorem  is  proved* 

We  now  obtain  a  simple  formula  for  the  sum   of  the   c. . 
This  quantity  arises  in  problems  of  aerodynamics  [2]  and  Its 
simplicity  allows  it  also  to  be  used  as  a  check  on  some 
alleged  inverse. 

Corollary  „ 


(2.9) 

•J 


TZ  °ij-=  XZ  (&k-bk)  = 


Proof o   If  we  apply  (2<,5)  to  p(x)  =  B.(x)   and  set 

J 

x  =  b,   we  obtain  that 


(2.10)         J__   B  (ai)Ai(bk)  =  6.k   . 

By   symmetry   this    is    also   valid    if      A,a      and      B,b      are    inter- 
changed.     Thus   we   get   that 

s   =     r~  a,.   CA.(b.)B.(a.)    -  J~  b.    YZ  A  j  (b.  )B±  (a  ^ 

=     IZI  a     -    HI  b.       . 

j         J  i 

For   the   special  case  of  the  Hilbert   segment    (1<>2),   formula 
(2«9)    gives   Smith's    [6]   formulas      s   =  n(p+n)      and   for   the 
matrix   of    (1„3)   we   get    the   formula   of  Collar    [2]s      s    =  pnc 

3 .      Row   and   column  sums. 

Let    the  row   and    column  sums    of      G-     be,   respectively, 

_   7   . 


(3.1)         SZOij  -«!.  i-c^  -Pj 

and  define  the  diagonal  matrices   ^  »D„   by 

Da  =  £ai>a2»°°->an] 
DB  =  [Pi»P2^°°'pn] 
Theorem  2„   The  matrix  >H   then  satisfies  the  following 
relations  I 


-A(b1)  B(a.) 

f— ^~  ,   P«  =  L 

B(b.)      J    A»(a.) 


(3o2)        a  = r-^—   ,   B.  = 


(3.3)  H"1  =  DaHTDp  . 

Proofs      Assuming    (3*2)    to  be   true  we  note    that      c.. 
can  be  written   in   the  form 

„   M  -         1  ^V  B(a.1}     =       ^ 

U'4'  °IJ     "bi-aj      B(b7"     ~^jT  aj-   bi 

which  gives    (3«3)° 

To  prove   the    formula  for      a.      in    (3<»2)    we  need   only- 
show   that 

B(a,) 

(3.5)  IZ    7- — r^TiTT"  =  1* 

j         (a.-  b±)   Av  (a..) 

However  for  any  function  f (x)   we  may  write 

,  _»      f(a.) 

(3.6)  f[a1,a2,..o,an]  =  6^"1Jf  (x)  =  ZZ     At(a.) 

where  the  left  side  of  (3.6)  is  the  (n-1)    divided 
difference  of   f (x)   with  respect  to  the  ak  (see  Milne- 
Thomson  [?]  pa  9)0   If  we  apply  (3*6)  to  the  polynomial 


f (x)  =  B(x)/(x-  b.)   of  degree   n-1  we  have  that 
6   "  f(x)  =  constant „   Since  the  coefficient  of   x 
in   f (x)   is  1,  (3°5)  is  proved,. 

An  alternate  proof  of  (3«5)  may  he  obtained,  without 
divided  differences,  by  noting  that  the  left  hand  side  of 
(3«5)  is  the  sum  of  the  residues  of  the  function 
B(x)/A(x) (x-  bi)   at  the   a . „   (This,  in  fact,  follows 
readily  from  (2.6).)   However  for  a  sufficiently  large 
circle   C   about  the  origin  in  the  complex  x-  plane 


2nii/n    (x- 


B(x) 


b.)A(x) 


dx  =  1 


which  proves  (3°5)° 

The  proof  for  the   B.   is  obtained  in  the  same  manner 
from  (3.5)'With  the  roles  of  A_,a   and   B,b  interchanged, 
and  the  theorem  is  proved  „   (The  formula  (3»5)>  incidentally, 
represents  an  extension  of  one  of  the  formulas  of  Collar  [2] 
for  the  sum  of  a  series.   Extensions  of  other  formulas  given 
in  [2]  can  likewise  be  obtained  from  (3<»6)  by  specializing 

fW.) 

!(..   Remarks » 

!•)   We  note  that  if  H   is  symmetric  then  a.  =  B. 


9  - 


and    if  we   set     D  =  D  ■   =  D„      then 

a  (3 


(k.l)  G     =     D  H  D 


2,)      Prom    (3»2)    and    (3.3)    one    immediately   obtains    the 
formula  for   the   determinant   of     H     up  to  a   sign.      That    is 

0|.2)  (det   H)2   =  Tl 


i       aiPi 

and  for  H  symmetric  we  get  that 

n(n-l)      A1  (a.) 

(U.3)     det  H  =   (-1)    2    Tl     ~   . 

i    B(&i) 

The  signs  may  readily  be  determined,  by  induction,  by  using 
(2.2)  and  the  formula  for  the   (n+l,n+l)   element  of  Gn+i: 

c  .-,       ,n  "det  H  .  n   =  det  H 
n+l,n+l      n+1  n 

where  Hn  =  H  and   Hn+1  ={a  -b  )   wlth  1^i'J^n+l0 

3o)   Formula  (2,l±)   may  also  be  applied  to  the  problem 

of  obtaining  a  least-squares  fit  of  a  function  f (x)   on, 

__    *± 
say,   (0,1)   by  a  function  of  the  form  >  y.x      .   (All 

i   x 
variables  are  here  assumed  to  be  real.)   The  normal  equations 

for  this  problem  yield  a  matrix  of  the  form  (1.1 )  with 

b.  =  -a.-l,  and  the  Y*    may  be  obtained  by  applying   G   to 

a . 
the  moments  of  f (x)   with  respect  to  the   x 
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For  fitting  a  function  of  two  variables  the  problem 
to  determine  y .  .   such  that 

11  —  a.      a.    « 

J   /    [f  (x,y)    -  ZL  Yli   x    x  y    xr   dx   dy   =  min 

0   0  i,j 


is   solved   by 


K     =     D  H  D  F  D1    H1  D1 


Where      H  "  ia1+a.j+l   3  •        H'    =   laUal-M      }>      K    =   ITij}- 

(l   1  a,      a'  ) 

P  =     J    J    f(x,y)    x       y   J    dx  dy      ,      and     D,D        are   the   diagonal 

lo  0  J 

matrices    corresponding   to     H,HV ,   respectively,    defined   above. 

Although  the  elements  c. .  get  quite  large  in  the  case 
of  the  Hilbert  matrix,  it  may  happen  that  for  suitable  choices 
of   the      a.      this   may  not   be   the    case. 

An  explicit   solution  for      y-      can  also  be    obtained   from 

(2.3)    for    the    equations 

1  a.        a' 

j(f(x)    -  JZL  Yj  x      )    x  K  dx    =  0,      k  =  1,2,. ,.,n 

0  *" 

ai 

given   the  moments    of      f (x)      with  respect    to   the      x         . 

\\a)      If  H   is  r.eal  and  symmetric  and  if  a.  >  b„, 
i  =  l,2,...,n  then  it  follows  from  (I.I4.)  that  all  the 
principle  minor  determinants  are  positive  and  H   is 
positive  definite.   In  this  case   a.  >  b.   for  all   I 


-  11  - 


ows 


and   j  5   that  is,  all  the  elements  of  H   are  positive, 
since  if  a.  <  a.,  then   0  <a.- b.  <  a  .-b.  =  a.-b..   Thus 
B(a.)>0  and,  since  A(x)  has   n  simple  zeros  at  the 
a.,  A'  (a.)   alternate  in  signD   Prom  (2.!j.)  it  then  foil 
that  if  a.  <  a .   for   i  <  i   then   (-1)1  ^  c. .  >  0  so 
that  then  G  has  the  same  checkerboard  distribution  of 
signs  as  the  inverse  of  the  Hilbert  segment.   From  (3*2) 
it  follows  that, in  this  case, the   a.   also  alternate  in 
sign  and  that   (-1)''   a.  >  0. 

J 

Let  \     and   u-  denote  the  smallest  and  largest  eigen- 

a,  -b. 

value  of  H  respectively,  and  let   M.  =  min   ( K 

k,k/^i        I  ak-  a±\ 

Then   it   follows   readily   that       |A.(b.)|    >  M1?"1  >  1      and    that 

1                        i 
X  <  min  <  min    ( r— ) ,    and    if   the      a.      increase   with     i 

"     i      cii  "     i        ai"bi  X 

then     min    (J^      (-1)1+J*    c.J"1   <   A.      We  get   for    the    P-condition 

[8]    of     H, 

^  >  max  (— -r— )  °  max  c.  .  >  max  [A.  (b.  )]  >  max  M7 
X  ~  i   ai_bi     i   il_i   Ti   "  i    1 

so  that  the   P-condition  of  H  may  get  very  large.   This 
number  has  been  estimated  for  the  Hilbert  segment  (102) 
with  p  =  0  by  Todd  [8], 

5o )   Prom  the  foregoing  we  may  also  get  a  simple 
proof  of  a  lemma  of  Ky  Fan  and  Gordon  Pall  [9]  which, 
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stated  in  a  slightly  more  general  form,  says  the  following; 
Let   a-,  ,a-> ,  o  o  o  ,a  ,,  ,  b,  , » .  <.  ,b  ,  t   be  colinear  complex 


^,*2 


n+l>  "1 


n: 


numbers  such  that  for  a  fixed   9   and   t   and  real  numbers 


py*i 


t   +Pje 


19 


b..    =    t    +  t,eie 


1   <    j    <  n+1 
1   <    J    <   n 


where 


ik'k) 


pl   >  ^1  >   p2  >  ^2  >  °  ° "    >  ^n  >    p 


n+1      . 


Then   there    exists    a   unique    solution   vector      x      to   the 
system 


(U.5)  Z 


n  x 


n 


+     x„ , ,    =  a .      s      1   <    i   <  n+1 
a.-  bn,  n+1  J      *  —   <j  — 


k=l      "J"    "k 


such   that 


2i9 


To  prove   this   we    let      u,    x",    a,   r      be   column  vectors    of 
order     n     such   that 

1 


uj  ~  lf   xj  "  xy   aj  *  ar   rj  :  an+1-  bj 

and  write   (l|.»5)  in  the  form 
(I4..6) 


1  <  J  <  n> 


T— 
Vr  x  +  x 


+1   :  an+l 


where   H   is  given  by  (lol)o   By  eliminating   xn+-i   we  obtain 
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T  —   — 
(U -7 )  (H  -  ur  )  x  =  a  -  uan+1   . 

A  single  computation  shows  that 

H  -  ur    =   LXH  L£x 

where   L,   and   Lp   are  diagonal  matrices  given  by 

Ll  =  [an+l-V"'e*an+l=an] 
L2  =  ^an+l  ~\'  •••»an+l'bn^ 

Thus  it  follows  from  (1+.7)  that 

x  =  L2H'1L^1(a-uan+1)  =  -LgBT^i  =  -L2a 

where      a      is    the   column  vector  with  components      a.      given 

by    (3.2). 

n  n 

If   we    let      R(x)    =     ~[\     (x    -    p.),    T(x)    =  TT    (x    -  tJ 

j=l  J  j=l  J 

and   let     R.(x),    T.(x)      be    the   corresponding  Lagrangian 

J  J 

coefficients    then 

A(b.)  RCt",)         ie 

J  B*  (bj)  OfCtj) 

This    implies    that   for      1   <   j    <  n 

Prom    (l+.lj.)    it    follows    readily    that   for      1   <    j    <  n 


R(f.)    <   o   <  t'(T,)      so   that      Y-    >   o 

J  J  <J 


11+    - 


In  the  particular  case  where   a.   and   b.   are  real 
and   t   and   6  are  chosen  to  be  zero  we  get  the  result 
as  stated  In  [9]  that  there  exists  a  unique   x  such  that 
x.>o  for   1  <  j  <  n. 
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